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The dielectric optical waveguide described in this paper has an annular 
cross section, the refractive index of which is higher than the indices of 
the material inside and outside of the ring. The solution of the eigenvalue 
problem is an approximation which is valid for small refractive index 
differences of the three media. The resulting approximate eigenvalue 
equation is far simpler than its exact counterpart. The cutoff conditions 
of the first three modes and the eigenvalue of the lowest-order mode are 
presented graphically. 

I. INTRODUCTION 

Dielectric optical waveguides have become interesting since it has 
been demonstrated that the losses of dielectric materials can be quite 
low. 1-3 The conventional optical fiber waveguide consists of a solid 
dielectric core surrounded by a dielectric material with lower refractive 
index. This structure has the advantage of being particularly simple. 
In this paper we analyze a different type of dielectric optical waveguide. 
Our structure consists of a dielectric tube which is filled and surrounded 
by dielectric material with lower index of refraction. Such structures 
have been analyzed before. 4 A complete discussion of and references to 
the literature can be found in Ref. 5. However, the exact analytical 
treatment is extremely complicated so that the resulting theories are 
hard to apply to practical situations. Recently A. W. Snyder 6 has 
shown that the analysis of the conventional optical fiber waveguide 
becomes much simpler if one assumes that the difference of the refrac- 
tive indices of core and cladding material is only slight. D. Gloge 7 
has developed Snyder's theoretical work even further. On the basis of 
this simplified method of analysis, it is possible to treat rather com- 
plicated optical waveguide structures more simply and still obtain 
results that are applicable to most cases of practical interest. The 
simplification that results from the assumption of only slight index 
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differences is particularly appropriate since most practical dielectric 
optical waveguides satisfy this requirement. 

An approximate eigenvalue equation for the modes of the tube 
waveguide and an equation describing the cutoff conditions for these 
modes are derived in this paper. The connection of the tube waveguide 
with the dielectric slab waveguide in the limit of large tube radii and 
with the solid-core optical fiber in the limit of small tube radii is 
discussed. The cutoff condition for the first three modes is represented 
graphically and the eigenvalue of the first mode is plotted. 

For simplicity we assume throughout that the index of the material 
inside of the tube is smaller than or equal to the index of the material 
surrounding the tube. 



II. APPROXIMATE MODE ANALYSIS OF THE TUBE WAVEGUIDE 

The geometry of the tube waveguide is shown in Fig. 1. Following 
the ideas of Snyder and Gloge 6 - 7 we express the electric and magnetic 
field components in rectangular cartesian coordinates and write 8 



E x = 



dE z dH.> 

+ "M 

dx dy t 



(1) 



B v = 
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dE 2 

dy 



OJ/i 



dH t 

dx 



(2) 




Fig. 1 — Cross section of the tube waveguide. 
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It is assumed that the time and z dependence of the field are determined 
by the factor 

««••-*>. (5) 

/3 is the propagation constant in z direction, w is the radian frequency. 
The parameter Kj is defined as 

K) = n % 2 - P 2 j = 1, 2, or 3 (6) 

with 

k> = »« fW i # . (7) 

Equations (1) through (4) are exact and are simply four of Maxwell's 
six equations written in a different form. If the refractive indices are 
all nearly the same, ni « n 2 « n 3 « n, we must have 

- nfc, (S) 

so that 

Kj « 0. (9) 

It is thus apparent that the transverse field components are much 
larger than the longitudinal components if the refractive indices are 
all very nearly the same. The longitudinal field components are 
obtained as solutions of the wave equations. 8 We use the expressions 

E z = [Z, + 1 (K } r) sin (u + 1)0 + Z,_i(A» sin (v - 1)0} (10) 

2n,jk 

VT iAjKj 
—rr- \Z, + i(K 3 r) cos („ + 1)0 
Ho 2k 

- Z^ x {K,r) cos (v - 1)0} . (11) 

The exponential factor (5) has been suppressed. Z v {Kjr) is a cylinder 
function and A,- is an arbitrary amplitude factor. Substitution of (10) 
and (11) into (1) through (4) results in 

E y = AjZ y (Kjr) cos u<f>, (12) 

[to 

H x = - UjAjJ— Z v (Kjr) cos v<t>, (13) 

Vmo 
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and 

E x = 0, H y = 0. (14) 

These equations are approximations that hold provided that (8) is 
applicable. 

In the three regions shown in Fig. 1, we use the following cylinder 
functions and parameters, assuming ni > n 2 > n 3 , 



-iK a = 6 = (/3 2 - n,fc»)*l 
Z, = J v (i9r) I 


for 


^ r ^ 6 


(15) 


K x = k = (n\k* - /S 2 )*l 

Z y = J v ( K r) + BN v (kt) 1 


for 


6 ^ r ^ a 


(16) 


-iK 2 = 7 = (P 2 -n 2 V)*J 
Z, = Hr tt) (tyr) 1 


for 


a ^ r ^ °o . 


(17) 



J v , N v are the Bessel and Neumann functions and H r (1) is the Hankel 
function of the first kind. We use the functions with imaginary argu- 




Fig. 2— The cutoff value of the normalized frequency V for G = as a function 
of the ratio of wall thickness d to tube (outer) radius a. The cutoff value for the 
mode v = is V = 0. 
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Fig. 3 — The cutoff value of V as a function of d/a for G = 0.7. 



ments instead of modified functions. This practice is in agreement with 
Jahnke and Emde 9 as well as with Gradshteyn and Ryzhik. 10 The 
amplitude coefficient Aj appearing in (10) through (13) also assumes 
different values in the three regions of space, j = 1, 2, 3. 

In order to be able to match the boundary conditions, we transform 
the cartesian components to vector components in cylindrical polar 
coordinates 



Ef = — E x sin <f> + E u cos 

= ^AjZ y (Kjr)lcos (p + 1)0 -f- cos (v - 1)0] 



(18) 



H* = \njA ix \— Z,CK»[sin {v + 1)0 - sin (v - 1)0]. (19) 

The boundary conditions require that E z , E^ H z , and i/ are con- 
tinuous at r = b and at r — a. The boundary conditions thus provide 
us with eight equations. However, our approximate treatment of the 
modes yields only four arbitrary constants, A\, A 2 , A 3 , and B. The 
situation is further aggravated if we observe that the field components 
contain the functions sin(i> + 1)0 and sin (i» — 1)0 and corresponding 
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expressions with the cosine function. Since the boundary conditions 
must hold for all values of <f>, we must equate the coefficients of the 
functions with v + 1 independently of those with v — 1, thus doubling 
the number of equations. In spite of this apparently hopeless situation, 
approximate solutions are possible. The equations that result from 
the requirement that the tangential components of E are continuous at 
the boundaries differ from the corresponding equations resulting from 
the boundary conditions for H only by factors ny. Since our approach 
is based on assuming that all n y are nearly the same, we set these 
factors all equal to the same average value n and, after dividing by n, 
obtain equations that duplicate those obtained from the boundary 
conditions for E. To the approximation considered here, the E and H 
boundary conditions lead to the same set of equations, reducing their 
total number to one half of the original number. To satisfy the con- 
tinuity of E z and E* at r = b, we equate the coefficients of sin (v + \)<t> 
and obtain 

J, (Kb) + BN,(Kb) 

A 3 = — Ai (20) 

J,(i0b) 

and the determinantal condition 

J, +1 (idb) J,+i(Kb) + BN y+1 (Kb) 

id = K- — — \.*1) 

J,(i6b) J V W + BN,(Kb) 

Equating the coefficients of sin (v - 1)0, we obtain a similar result 
with the only difference that v + 1 ifl replaced by v — 1 in (21). 
Equation (20) remains the same. With the help of the recursion 
relations 9 ' 10 for cylinder functions, it is easy to prove that eq. (21) is 
not changed if v + 1 is replaced by v - 1. We thus find ourselves in 
the happy position of being able to satisfy-at least approximately-the 
eight boundary conditions at r = b with the two available coefficients. 
In an analogous fashion, we obtain from the boundary conditions at 
r = a the equation 

J^ + OfM 
HW(na) 

and the determinantal condition 

. Hl+i(iya) _ ,7,,+iM + ■BiV y+ i(/ca) 
17 HW(iya) ~ * /,(<ca) + BN v (k(i) 
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Fig. 4— Same as Fig. 3, G = 1.5. 

The eigenvalue equation is obtained by eliminating B from (21) and 
(23). We find 

Kj y (ieb)J v+ i(Kb) — i6J v (Kb)J, + l (idb) 
Kj y (idb)N,+i(Kb) - idN ',(«&)./ \+i(idb) 



kH v (iya)Jy+i(Ka) — iyJ v (Ka)H>+i(iya) 
" KH^(iya)N v+1 ( K a) - iyN v {Kd)H^(iya) 

The constant B is given by 

Kj,{i6b)J, + l (Kb) - iBJ,(Kb)J v+l (idb) 



B 



(24) 



(25) 



kJ \(i6b)N ;+i(ic6) - iON>(Kb)J, +1 (i6b) 

It can be shown that the eigenvalue equation (24) specializes to the 
eigenvalue equation (d = a — b) 

k(7 + 0) 



tan nd = 



yd 



(26) 



of the asymmetric slab waveguide in the limit a — > <» and 6 — > °o . The 
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distinction between TE and TM modes of the slab waveguide is lost 
in our approximation. 

We are mainly interested in the cutoff conditions of the modes in 
order to find the separation between the lowest and the next higher 
mode which determines the range of single-mode operation. Cutoff is 
defined by the condition 

7=0. (27) 

In order to determine the cutoff frequency, we use (27) in (24) and 
obtain the cutoff equation 



■'■{ iGv -yA v ^)- iGN '( vb i) j A iGV l 



(28) 



The parameters appearing in this equation are defined as follows: 

V = (n! - nlykd (29) 



and 

o = /n 



2 2 vi 

— T),\ I 



—]■ (30) 



For future reference, we derive an interesting relation for the mode 
with v = 0. The cutoff condition (28) allows the solution V = for 
v = 0. We find as the condition that V = 0, the relation 



2 2 



1 3 

In the limit of infinite radii, a — > °o and 6 — » « , we obtain from (26) 
the cutoff condition for ni > n 2 ^ n 3 

7 = arctan G. (32) 
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Fig. 5— Same as Fig. 3, G = 5.3. 



III. DISCUSSION AND NUMERICAL RESULTS 

We restrict our discussion to the case n 2 ^ n s . The cutoff values of 
V [see eq. (29)] as functions of d/a are plotted in Figs. 2 through 6 
as solid lines for the first three modes, i» = 0, i»=l, »» = 2. The V 
parameter characterizes the tube waveguide in the limit of large radii 
where it can be regarded as an asymmetric slab waveguide that is 
bent into a circle. For large values of d/a, the waveguide approaches 
an optical fiber with solid cylindrical core. Such a fiber is characterized 
by a different V parameter which can be expressed in our notation as 
(a/d) V. In the limit d/a = 1, both parameters coincide. The parameter 
(a/d)V is plotted in the figures as a dotted line only for the mode 
v = \. Each curve, Figs. 2 through 6, is plotted for a different value 
of G [see eq. (30)]. The relation between the G values and the ratios 
of ni/tij and ni/n 3 is shown in Fig. 7. 

Since the tube waveguide approaches the asymmetric slab in the 
limit d/a = and the solid-core fiber in the limit d/a = 1, the cutoff 
values V = V e can be predicted at these limits. For d/a = 1, the 
cutoff value of the solid-core fiber is obtained from 



J,-i{V e ) = 0. 



(33) 
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The values of V c at d/a = 1 are, therefore, V c = for v = 0, V c = 2.405 
for v = 1, and V c = 3.83 for v = 2. 

At d/a = 0, we find the cutoff values for V from the asymmetric 
slab formula (32). 

We see from Figs. 3 through 6 that the curve with v = seems to 
end on the horizontal axis. Actually, the curve must be continued along 
the horizontal axis at V = to the point d/a = 1. However, to the 
left of the point where the v = curve reaches the horizontal axis, 
V = is not a legitimate solution of the cutoff equation. The d/a 
value at which the v = curve touches the horizontal axis is given by 
(31). This phenomenon finds the following simple explanation (private 
communication from E. A. J. Marcatili). As the mode approaches its 
cutoff value, we have 7 = 0, and also 6 = 0, and 7=0. There is no 
longer any transverse field variation so that the field "senses" only 
the average value n 2 of the dielectric constant of the tube, 

* = ^' + ^'-"^'. (34 ) 

7ra 2 

The mode is no longer guided if the average dielectric constant n 2 is 
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Fig. 6— Same as Fig. 3, G = 9.6. 
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Fig. 7 — The parameter G is shown as a function of (/ii/n 3 — 1) for several values 
of ni/nj. 



equal to the dielectric constants n\ of the outer medium. The condition 
n 2 = n\ also leads to (31). 

Figures 2 through 6 show that the cutoff points for all three modes 
(v = 0, 1, and 2) approach each other arbitrarily closely for small 
values of d/a. It is thus apparent that single-mode operation in the 
mode v = becomes increasingly more difficult to achieve as the 
ratio d/a is decreased. 

Finally, Fig. 8 shows a plot of the transverse decay parameter ya 
as a function of d/a. The curve is intended to provide information 
about single-mode operation of the tube waveguide. However, since 
it is desirable to operate with a tightly guided mode field, we chose as 
the operating point for V the cutoff value of the second mode, V = V ex- 
it is a remarkable feature of the curves of Fig. 8 that they become 
independent of d/a for small values of this ratio. However, from the 
point of view of field confinement, we must remember that it is not 
ya but rather yd that characterizes the rate at which the field decays 
from the guiding structurc-the tube-in radial direction. The crosses on 
the curves in Fig. 8 indicate the points at which we have yd = 1, the 
circles indicate the points yd = 0.5. 



434 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1973 



1.2 



Ya 



1.0 



0.8 



0.001 



k = 

v = v c , 

X Vd= 1 
O yd = 0.5 



I I I I 




J I L 



J L 



0.01 



1 



d/a 



Fig. 8— This figure shows the parameter ya as a function of d/a for several values 
of G. V is taken equal to the cutoff value V cX of the mode v = 1. 



IV. CONCLUSIONS 

The eigenvalue equation for the modes of the tube waveguide has 
been derived with the help of an approximate technique. The cutoff 
values of the normalized frequency parameter V [denned by (29)] are 
presented for the first three modes of the tube waveguide and the 
transverse decay parameter ya is plotted for the lowest-order mode 
under the assumption that this mode is operated at the point of cutoff 
of the next-higher-order mode. 

One may wonder if it is possible to use the tube waveguide with 
single-mode operation to increase the mode radius compared to the 
mode radius of the lowest-order HE U mode of the conventional fiber. 
This question has been studied under the assumption that both wave- 
guides provide an equal amount of field confinement. The criterion for 
field confinement is the radial decay parameter y. It was consequently 
assumed that the y values of both waveguides are identical. Surpris- 
ingly, it was found that under the condition of single-mode operation 
for each guide and the requirement of equal field confinement, the 
mode radius of the tube waveguide is only slightly (approximately 
40 percent) larger than the radius of the HEn mode of the optical 
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fiber. A large mode radius appears desirable from the point of view 
of alignment tolerances of fiber splices. With a large mode radius, the 
alignment of the waveguide centers would be less critical. The tube 
waveguide does not appear to offer a significant advantage for this 
purpose. 
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